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Doubly fed Induction Generators (DFIGs) are quite common in wind energy 
conversion systems because of their variable speed nature and the lower 
rating of converters. Magnetic flux saturation in the DFIG significantly affect 
its behavior during transient conditions such as voltage sag, sudden change in 
input power and short circuit. The effect of including saturation in the DFIG 
modeling is significant in determining the transient performance of the 
generator after a disturbance. To include magnetic saturation in DFIG model, 
an accurate representation of the magnetization characteristics is inevitable. 
This paper presents a qualitative modeling for magnetization characteristics 
of doubly fed induction generator using neuro-fuzzy systems. Neuro-fuzzy 
systems with one hidden layer of Gaussian nodes are capable of 
approximating continuous functions with arbitrary precision. The results 
obtained are compared with magnetization characteristics obtained using 


Saturation discrete fourier transform, polynomial and exponential curve fitting. The 
error analysis is also done to show the effectiveness of the neuro fuzzy 
modeling of magnetizing characteristics. By neuro-fuzzy algorithm, fast 
learning convergence is observed and great performance in accuracy is 
achieved. 
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1. INTRODUCTION 

Wind energy is gaining importance in the present-day power scenario due to the rising concerns of 
the environmental issues along with the depletion of conventional sources and the increase in power demand. 
As reported by the world wind energy association (WWEA) that the total wind capacity in the world has 
attained 486,661 MW by the end of 2016. This report also reveals that 54,846 MW were added in 2016. This 
indicates a growth rate of 11.8 % in 2015 [1]. The main part of the wind energy conversion system is 
electrical generator. Nowadays wind turbines widely use doubly fed induction generator (DFIG), as it can 
operate for variable wind speeds. The stator of doubly fed induction generator (DFIG) is straightly connected 
to grid whereas the rotor is connected to the grid through a bidirectional voltage source converter. 
Accordingly, the control of the machine can be implemented with the converter of power rating of 25-30% of 
the rated power of the machine. This results in less cost, size and weight of converter and also less losses and 
increase in efficiency [2]. 

Wind energy conversion systems based on doubly fed induction generator are highly susceptible to 
grid disturbances such as voltage dips. Earlier DFIGs were disconnected from the gird during voltage dips. 
Due to the increased use of wind power, this disconnection of wind turbines may lead to loss of power and 


Journal homepage: http:/iaescore.com/journals/index.php/IJECE 


24 o ISSN: 2088-8708 


also to the system stability problems [3]. Besides, the sudden drop in grid voltage may cause over voltages 
and over currents in the rotor due to the magnetic coupling between stator and rotor. Due to these high 
currents, thermal breakdown of converter may occur. So according to new grid code requirements, all wind 
turbines must remain connected in the grid during short circuit also [4]. 

When the machine is subjected to several transients such as loading and short circuit faults, the 
ferromagnetic core of the machine saturates resulting in significant change in the electromagnetic parameters 
of the machine. Hence, actual machine responses can be different from the ones predicted by machine models 
that assumes linear magnetization characteristics. Therefore, the inductances that are usually treated constant 
in the conventional model are functionally dependent on machine currents. To have a machine model that 
resembles closely a real doubly fed induction generator, new terms representing the main flux saturation have 
to be included. Eventhough several works have done in the field of saturated induction motors, the effect of 
the magnetic saturation is usually neglected for DFIG. The saturation effect in the induction motor is 
analytically modeled in either the stator and rotor leakage flux paths [5]-[6] or in the main flux path [7]. 
On the other hand, experimental methods are utilized in [8] to incorporate the magnetic saturation in the main 
flux path, and in [9]-[10] to compensate for both the main and leakage flux paths. In [11] the main flux 
saturation is developed by replacing the unsaturated magnetizing inductance or reactance by its 
corresponding saturation function. Moreover, recent works are also available in [12]-[14] which considers the 
saturation based on analytical methods. 

The effect of magnetic saturation in an electrical machine is realized by the mathematical 
formulation of the magnetization characteristics. Hence, an algorithm that can show functional relationship 
between current and flux will be useful for including the magnetization phenomena in electrical machine 
model [15]. Various methods have been presented by researchers to model the magnetizing characteristics 
from experimentally determined data points. It is found that functions such as exponential, rational and 
polynomial curve functions are used in literature to fit the magnetizing characteristic data points [16]-[24]. 
Soft computing techniques such as feed forward multilayer ANN networks are also developed to model 
saturation in electrical machines [25]-[27]. 

The objective of this paper is to develop a neuro-fuzzy model to represent the magnetization 
characteristics of the machine. As a universal approximator, neuro-fuzzy systems with Gaussian functions is 
used for the function approximation of magnetization characteristics. The mean square error is calculated to 
compare the results with the frequently used regression models such as polynomial, exponential and discrete 
Fourier transform (DFT). Neuro fuzzy algorithm is advantageous as its recursive nature of the computation 
results in enhanced function accuracy and also the parameters needed for approximation are only the centre 
and the width of the Gaussian membership function. Hence programmability for parametric learning is easy 
and straightforward. Even though, the polynomial models are advantangeously linear in parameters but it is 
limited in terms of generalization capabilities. As DFT is a sinusoidal approximation method, to increase the 
accuracy, order of cosine terms should be increased. This work can be considered as a preliminary step to a 
future implementation of neuro-fuzzy approximation in the modeling of DFIG for transient stability studies. 

The section 2 explains the model of doubly fed induction generator and section 3 describes the 
nonlinearity of magnetizing characteristics and section 4 three existing models used to validate the proposed 
model and section 5 describes the general neuro-fuzzy structure and application of the neuro-fuzzy algorithm 
to develop the magnetization model of DFIG. The section 6 describes the numerical validation using the 
above algorithms. The section 7 includes the results and discussions. 


2. DOUBLY FED INDUCTION GENERATOR MODEL 

Doubly fed induction generator is essentially a wound rotor induction machine with its stator 
directly connected to the three-phase grid and rotor to the grid by two back to back converters namely the 
rotor-side converter and the grid-side converter [28]. Schematic diagram of DFIG as shows in Figure 1. 
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Figure 1. Schematic diagram of DFIG 


The DFIG is modeled in synchronous dq reference frame for the reason of simulation, analysis and 
control. Assuming the motor convention and taking into account of the flux linkage transients, the voltage 
and flux equations of DFIG are represented as [28]. 


P d 
Vsq = Rsisq a an Psa + OsWsq 
r? d 
Vsa = Rsisa + g Usa = WsWsq 
š d 
Urg = Rrirq + q Pra + (ws — wr)Pra 


i d 
Vsq = Rsisq + an sa — (ws — 0, )Wsa (1) 
The stator and rotor fluxes are given by 


Wsq = Lsisq + Lmirq 
Psa = Lsisa + Lintra 
Pra = Lrirq + Lintsq 
Pra = Lsira + Lmisa (2) 


Where Ls = Lis + Lm 
Ly = Ly + Ln 


3. NON-LINEARITY OF MAGNETIZATION CHARACTERISTICS 

The saturation of the ferromagnetic material used in the construction of electrical machines is 
responsible for the nonlinearity of magnetization characteristics. The conventional model of doubly fed 
induction machine shown in (1)-(2) employs constant value of magnetizing inductance as it assumes linear 
magnetization characteristics. As the flux is nonlinearly related to the magnetizing current, the conventional 
DFIG model formulated on the basis of linear magnetization characteristics fails to obtain precise results and 
causes performance predictions nearly inaccuratein many instances. Hence, more precise study of 
magnetization characteristics is necessarily needed to include the effect of magnetic saturation in the machine 
modeling [29]. Figure 2 shows the non-linear relationship between the magnetizing current and the 
opencircuit voltage. The main flux which is a function of magnetizing current can be expressed as 


Ly = Yrm) 3) 


lm 


As far as steady state is considered, the saturated magnetizing inductance is proportional to the chord slope of 
the saturation characteristic. 
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In transient analysis of DFIG the incremental inductances are also required. 


dt di dt 
d di 
Ym = Lmt PA (4) 


The transient magnetising inductance Lmt is proportional to the tangential slope of the magnetization 
characteristic. The inductance Lm and Lmt varies with the magnetizing current and appropriate method to 
identify them is to express the interpolating curve of the magnetization characteristic as non-linear function. 
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Figure 2. Magnetization characteristic 


The dominance of magnetizing current depends on the load power factor. Hence, whenever 
magnetizing component of current is more, the core can have a tendency to saturate and the flux is no longer 
linearly proportional to magnetizing current. Then the degree of non-linearity of the magnetization 
characteristics has to be tracked. Under such condition, prediction about the success of operation is not 
possible without using numerical methods. Here lies the role of neuro-fuzzy method. Neurofuzzy system with 
back-propagation learning algorithms have the ability to learn a continuous nonlinear function from input and 
output data samples and effective in identifying nonlinear dynamic systems with unknown characteristics. 


Accounting for magnetizing flux saturation, voltage equations of DFIG can be rewritten as [29]. 


j ü dad. i å 
Vsq = Rsisq + Lis utsa + Lma a ma + Ws Listsa ay Linstma) 

3 d., d., ‘i š 
Vsa = Rstsa + Lis a lsd + Ling qma + ws(Lis İsa + Listing) 

r d d š P 
Vra = Rrirq + Lir a a + Lma a ima + (ws — @,)(Lisisa + Lmsima) 


r E d š 5 
Vra = Rrira + Ly ura T Lma dt lma — (ws a Wr) (Lining + Lms inq) (5) 


4. DIFFERENT METHODS TO FIT MAGANETIZATION CHARACTERISTICS 

To fit a given set of data points, different mathematical techniques are available in literature. Here, 
we propose neuro fuzzy algorithm for modeling the magnetization characteristics. The results obtained from 
neuro fuzzy method have been validated with those obtained by discrete Fourier transform, polynomial and 
exponential curve fitting methods. 
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4.1. Discrete Fourier transform method 

Discrete Fourier transformation estimation can be used for n measured data points Ym &Im to 
represent the magnetizing flux in sinusoidal format. The DFT of the magnetizing flux can be 
expressed as [30]. 


WUm) = ao + pea ajcos(wjlm) (6) 


where the coefficients are defined as 


ao = — f" di (7) 


Imax 
di 1 
= ag URE? [vie = Ik-1) — Vale kk-1)?] 


and 


2 


aj = — O) cos(wji)di 


EF 
ggk. sin(wjIk-1) . (Ik — Ik-1) + wp (608 (@ylx) — cos(wjlę-1)) + 
7 (sin(w;Ix) — sin(wjlIk-1))] (8) 


where Yg, wj and Imax are defined as 


Imax = Imn (9) 
Yk = We-Wr-1 (11) 
Ik-Ik-1 


The number of cosine terms defines the order of the DFT function. For obtaining more precise results, the 
order of DFT should be increased. 


4.2. Polynomial curve fitting method 
Polynomial regression is one of the most widely used curve fitting method. The function 
representing the magnetization data is expressed in the form of a polynomial of order j 


W(Im) = ao + Alm + a212 + + alj, = ao + Xl arl% (12) 


The coefficients that best fits the curve to the data points is obtained by error-least squares approach. 
The general expression for error using least squares approach is given by 


; 2 
Error = X} (Yi — (ao + Zl, ax1%)) (13) 
The coefficients of the best fit curve are the one that minimizes the error function. 


4.3. Exponential curve fitting method 
The most widely accepted curve fitting function is of the form given by 


WUm) = D2, a; ePm (14) 


The error function is given by 


e = D4 (W(Imj) — wi) (15) 
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By the method of least square errors, the best fit coefficients can be obtained by minimizing the errors 


Oe _ ðe _ ðe 0 
Oa, = Oaz ~ dan ~ 
de _ ðe ðe 
ðb ðb, — bp 


5. NEURO-FUZZY ALGORITHM 

Neural networks and fuzzy logic are two promising modern technologies. Neural networks are often 
regarded as “black box” approach since we cannot conceptually understand the pattern learned by the neural 
networks. Unlike neural networks, fuzzy system uses comprehendible linguistic terms and if-then rules [3 1]. 

The combination of neural networks and fuzzy logic led to the formation of neuro fuzzy system. 
A neuro-fuzzy system describes a fuzzy rule-based model using NN-like structure [32]. A collection of input 
output data is taken and a learning method comparable to the neural networks is employed to train and adjust 
the parameters of membership functions. A simple structure of ANFIS is presented in Figure 3 [33]. 


j ' 0 = 
Input | Rules | Output | ao utput 
Input | Membership | Membership = 
| Function | Function 


Figure 3. Structure of neuro-fuzzy system 


Nonlinear curve fitting with neuro-fuzzy algorithm can be used to obtain the relation between 
magnetizing current and flux linkages with precision. The block diagram of the method is shown in Figure 4. 
The error is the difference between the flux linkages obtained from the experiment and the interpolated ones 


Nonlinear Curve fitting 


i 
by neuro-fuzzy algorithm e 


Flux mapping 
interpolation 


Figure 4. Block diagram to obtain flux-current relation 
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5.1. Neuro-fuzzy algorithm to fit magnetization characteristics 

A neuro fuzzy model to fit the magnetization characteristics data points of DFIG is demonstrated in 
this section. The neuro-fuzzy system used here is a first order Sugeno System with single input and three 
guassian distribution membership functions. Three if-then rules are used. The back-propagation algorithm is 
used to train this network.The training is performed by taking magnetizing current as the input and open 
circuit voltage corresponding to the magnetizing flux as the output. The linear membership functions are 
chosen for output. 

The magnetizing current is the input which is represented by gaussian membership function. 
The linguistic variables for the fuzzy partitions of the magnetizing current are small, medium and large. 
The input magnetizing current is fuzzified using gaussian membership function is represented as [32]. 


mal) = expl- (57) (16) 


where x is magnetizing current and i stands for linguistic variables. The ancedent parameters are centers 
Ms Mm, Mı and variances Os, Om, O, of small, medium and large input membership functions. 

The back-propagation algorithm used to train neuro fuzzy model consists of two passes, forward 
pass and backward pass. In the forward pass, for a given input the actual output of the system is computed 
using following equation [36]. 


Yk = Diea vi (k)ci (17) 


where 
ay 


oj 


sso) 


vi(k) = (18) 


which represents the normalized firing strength of the rules and c; is the consequent parameter. In the forward 
pass, the ancedent parameters and the consequents that corresponds to the weights in the neural network are 
fixed. While during the backward pass, the error signal due to the difference of the real and the desired output 
of the model is propagated in backward direction and the ancedent and consequent parameters are updated by 
employing the gradient search algorithm. The error function at the kth iteration is denoted as [36]. 


J(k) = O, - 7.) (19) 


The error correction rules for c;,0;,,m; are given by 


cilk) = i(k ZE losci (20) 
ye ee (21) 
mm, (K) = mik- I mimio (22) 
where n1, N2, N3 are the learning rate parameters. 
BO) = URKO = —2e(k)vj(k) (23) 
UO — UE Sot WO — —2e(kyi (Lei = Diba via |S" (24) 
UO — $10 Set RO — _2eckyvi (kyle — Zit na] [=| (25) 


The Batch Mode Backpropagation Algorithm applied to Takagi Sugeno Kang (TSK) models with Gaussian 
Membership Functions is as follows [36]. 
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Step1: Set c; (0), m;(0) and o;(0) to small random numbers and set n=1. 

Step2: Compute v;(k) using equation for k=1,2,........... ,N, where m; and g; are replaced by m;(n — 1) and 
o;(n — 1) respectively; Compute f(k) using equation for k=1,2,........ ,N where c; is replaced by c;(n-1) 
Step3: From error signal e(k)=y(k)-9(k) for k=1, 2......,N; update c; , m; and o; using equations 
respectively; 

Step4: If some stopping criterion is satisfied, stop, otherwise set n=n+1 and go to step 2. 


6. NUMERICAL VALIDATION 

The ratings of doubly fed induction generator considered in this paper is shown in appendix [34]. 
The open circuit voltage and magnetizing current needed for neuro fuzzy modeling of the magnetization 
characteristics of doubly fed induction generator is taken from [34]. The method to obtain the magnetization 
characteristics is by performing the open circuit test at synchronous speed. A three-phase controllable 
amplitude power source at the rated frequency has to be applied to the synchronously driven doubly fed 
induction generator. Vary the applied voltage to measure the terminal current. As the DFIG is operating at 
synchronous speed, the rotor current is zero and the magnetizing current is the terminal current itself [34]. 

The main flux saturation characteristics approximated using neuro-fuzzy model is shown along with 
data, DFT, polynomial and exponential in Figure 5. The parameters of neuro fuzzy system obtained are 
mentioned in Table 1. The training is performed using back propagation algorithm. Figure 6 shows the input 
membership functions of neuro fuzzy system after training process. To validate the accuracy of the neuro 
fuzzy method, the discrete, polynomial and exponential curve fitting method are also applied to 
magnetization data points. The calculated results are found to be matching with the measured ones. 
The accuracy of the proposed model is assessed by calculating mean square error. The mean square error can 
be determined by the given expression. 


)-xX;]2 
MSE = ye el (26) 
where X(Imi) are the calculated values by the saturated levels at different magnetizing currents, 
Imi» Im2,lm3, = Imn, Whereas X1,X2........Xn are the measured data points at these magnetizing currents. 
The results shown in Table 2 illustrate the effectiveness of neuro fuzzy method in fitting the magnetization 
characteristics of doubly fed induction generator. Figure 7 shows the mean square error obtained using neuro 


fuzzy method while training at each epoch. The number of training epochs is set to 50 with an error tolerance 
of .0005. 
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Figure 5. Calculated main flux saturation characteristics of the DFIG using different curve fitting method 
along with data 
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Figure 6. Input membership functions 
Table 1. Parameters of the Neuro-fuzzy System 
Input Membership Function Parameters Output Membership Function Parameter 
Centre m Spread o Polynomial Constant 
MFI 2.1184 4.6507 -1.3506 
MF2 4.6874 3.5622 0.6567 
MF3 9.9767 1.6832 -0.0150 


Table 2. Mean Square Error of Different Methods 
Curve Fitting Method Mean Square Error (MSE) Maximum Error 


Neuro Fuzzy 0.0005 0.0402 
DFT 0.0088 0.0593 
Polynomial 0.4452 0.9304 
Exponential 0.4444 0.9101 
0.25 T T T 
0.2 F J 
0.15 f J 
S 
wi 
oif 3 
0.05 F J 
0 | 
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Figure 7. Training error using Neuro-fuzzy system 
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7. CONCLUSION 

Neuro Fuzzy System for modeling the magnetization characteristics of doubly fed induction 
generator has been developed in this paper. The propriety of using this model in estimating the magnetization 
characteristics of doubly fed induction generator has been investigated. Using predetermined magnetizing 
flux—current data sets, a neuro-fuzzy system can tune the parameters of membership function by using back 
propagation algorithm. The non-linearity of the magnetization characteristics is effectively plotted by this 
method. The mean square error has been calculated to measure the accuracy of the proposed method. Results 
obtained are compared with the discrete Fourier transform, polynomial and exponential regression methods. 
This shows the effectiveness of the proposed method in the modeling of magnetizing characteristics. The 
neurofuzzy model is also easier to implement as it needs to tune only parameters of input and output 
membership function. 

It is concluded that neuro-fuzzy model developed can be used for magnetizing modeling of DFIG 
for its accuracy and reliability. Therefore, as the next step of the research, the magnetization model 
developed can be included in the state space model of DFIG. This will lead to a comprehensive and precise 
model which can be used for the performance analysis of doubly fed induction generator. In this sense, this 
work can be seen as a study previous to a future implementation of neuro-fuzzy approximation in the 
modeling of DFIG for steady state and transient stability studies. 
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APPENDIX 
DFIG Machine Ratings 
Machine Ratings 

Rated Power 2.78 KVA 
Rated Voltage 208 V 
Rated Frequency 60 Hz 
Rated Power Factor 0.72 
Rated Speed 1,720 rpm 
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